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Abstract 

We consider Freund-Rubin-type compactifications which are described by (p+g)-dimensional 
Einstein gravity with a positive cosmological constant and a g-form flux. Using perturbative 
expansions of Kinoshita's ansatz for warped dSp x S'^ and AdSp x S'^ spacetimes, we obtain 
analytical solutions describing the warped branches and their respective phase spaces. These 
equations are given by inhomogeneous Gegenbauer differential equations which can be solved 
by the Green's function method. The requirement that the Green's functions are regular 
provides constraints which determine the structure of the phase space of the warped branches. 
We apply the perturbation results to calculate the thermodynamic variables for the warped 
dSp X S'^ branch. In particular, the first law of thermodynamics can be reproduced using 
this method. 



^Email address: phylyk@nus.edu. sg 



1. Introduction 

Compactified solutions for de Sitter (dS) and anti-de Sitter (AdS) spacetimes have been 
the subject of interest from various viewpoints. For instance, (p + g)-dimensional spacetimes 
of product form AdSp x Aig are frequently considered in the context of supergravity and 
string/M-theory, as well as being an important system for the AdS/CFT correspondence. On 
the other hand, as scenarios of inflationary cosmology and the present day acceleration are 
approximately modeled by the de Sitter spacetime, then dSp x Aiq comes under important 
consideration when one attempts to incorporate cosmological models into higher- dimensional 
theories of (quantum) gravity. Thus it is of much interest to understand the physics required 
and caused by such scenarios. 

One of the important questions to address is the stability of such systems. For instance, 
the question of stability is important for cosmology as it is well-known that the effects of 
cosmological perturbations from matter fields play a significant role in the early universe. 
In the context of higher- dimensional supergravity, Freund and Rubin [1] have shown that 
the existence of a higher-dimensional flux will induce a natural compactification of the extra 
dimensions. If the spacetime is maximally symmetric, the compactified solution is of the 
form AdSp x S'^, where p and q are respectively the number of macroscopic and internal 
dimensions. This was explored further by De Wolfe et al. [2] who studied the stability of 
general AdSj, x Aig configurations, where A^^ is an Einstein space with positive curvature. 
In these theories, the negative curvature of the external spacetime is due to the presence of 
a g-form flux in the internal space. 

If a positive bulk cosmological constant is included in the theory, more interesting cases 
can arise. For a sufficiently small flux strength, the curvature of the external manifold can be 
positive and hence describes the de Sitter spacetime. Since present cosmological observations 
and early inflationary epochs are modeled by the de Sitter spacetime, this allows us to 
consider flux compactiflcations within cosmological models. For example, such a model has 
been considered by Carroll et al. [3] to construct a mechanism where the dSp x 5*^ geometries 
can dynamically arise from a spontaneous compactiflcation of an initially uncompactifled 
{p + q')-dimensional background. 

As mentioned above, a problem frequently considered in cosmology is that of pertur- 
bations in cosmological models, which are primarily studied to understand the formation of 
large scale structure in the universe and also to investigate the stability of compactiflcations 
in cosmological solutions. This was considered in the context of Freund-Rubin-type com- 
pactiflcations by Martin [1], which is related to earlier work by Bousso et al. [5], Contaldi 
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et al. [6] and Frolov [7]. Martin's result is that for small flux, conflgurations in the form 
dSp X 5''' are unstable against homogeneous perturbations, and beginning in the case (? = 4, 
for large flux, higher modes of perturbations (/ = 2, 3, . . .) begin to enter. For the case 5' > 5, 
there are no more stable modes for dSp x S'^. 

The onset of the threshold unstable mode may signal a presence of a new branch of 
solutions which are possibly non- uniform. In Ref. [8], Kinoshita looked for one of these 
solutions by generahzing dSp x S'^ into an ansatz corresponding to a warped product of the 
form 

ds^ = A{yfg^,{x)dx^dx'' + g^niy)dy"-'dy^ , (1.1) 

where x'^ denote coordinates of the external spacetime and y"^ denote coordinates in the 
compact space. Using this ansatz, a new branch was found numerically for q = 4 where 
the compactifled space is a deformed sphere. More recently, Minamitsuji and Uzawa [9, 10] 
obtained warped de Sitter compactiflcations with the topology dSp x x S'^~^. In contrast 
to the one studied by Kinoshita, exact solutions were obtained for this case. 

Interestingly, the results reviewed above parallels the recent developments regarding the 
stability of higher- dimensional black holes. In 2003, Emparan and Myers [11] showed that 
ultra-spinning Myers-Perry black holes are subject to an instability, which led to a proposed 
construction of a black hole phase diagram whereby a new branch of "lumpy" black holes 
are postulated to emanate from the Myers- Perry branch [12, 13]. 

Subsequently, numerical analysis on various cases of Myers- Perry black holes have iden- 
tifled the points of marginal stability (which are points where the new branches are expected 
to appear) in the respective phase diagrams of singly-spinning black holes [14, 15], and black 
holes with equal angular momenta [16]. Due to the extreme complexity of the Einstein 
equations for such systems, most of the results were obtained using advanced numerics and 
heuristic constructions. 

The Freund-Rubin equations we are considering here have a slight advantage over black 
holes due to the fact that the instability appears in the scalar sector with respect to dSp x S'^, 
thus making it more tractable to analysis. The introduction of warping in Ref. [8] singles 
out only one direction where S'^ symmetry is broken, while preserving the others. Such a 
system was solved numerically by Kinoshita to demonstrate the existence of a new branch 
of warped solutions for p = A, q = 4, and a phase diagram is obtained.^ Subsequent work 

^The Freund-Rubin phase diagrams considered here will be parametrized by the Hubble constant and 
flux strength. 
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by Kinoshita and Mukhoyama [18] shows that the entropy of the warped branch is higher 
than the unwarped branch whenever the latter suffers from instability, providing further 
support to the correlated stability conjecture by Gubser and Mitra [17] which was originally 
introduced for black branes with infinite translational symmetry. 

In this paper, using perturbative methods, we obtain analytical expressions to describe 
the phase diagrams for the warped branches. As we will show below, in the phase space, 
we are able to write down the equation for the warped branch as a series expansion about 
the point of marginal stability. The result we obtained supports the numerical analysis of 
[8], particularly, for g = 4. Continuing the analysis to g = 5, 6 and 7 we calculate the full 
phase diagrams numerically and also find agreement with the perturbative analysis. This 
procedure is analogous to the work by Gubser [19] in the study of non-uniform black branes. 
However for black branes even the perturbative expansions require numerical solutions. In 
our case, the equations are sufficiently simple and analytical solutions can be obtained. By 
analytical continuation of the Hubble parameter, the perturbative analysis can be easily 
extended to include anti-de Sitter spacetimes. Thus both dSp x S'^ and AdSp x S'^ solutions 
can be treated equally well."^ 

This paper is organized as follows: in Section 2 we review the theory for warped com- 
pactifications with flux, as well as setting up and defining the notations to be used throughout 
the rest of the paper. In the same section, we will also reiterate the main results of Refs. 
[5, 4] regarding the stability of dSp x JUg. In Section 3, we present the perturbative analysis 
of Kinoshita's ansatz. In particular, we will focus on the even modes (/ = 2N where is an 
integer) of the perturbative solutions. We show that at first order, they are equivalent to the 
results in [5, 1]. Proceeding to second and higher orders we obtain the analytical structure 
of the phase diagram near the points where the warped branch emanates from the unwarped 
branch. There is a subtlety in the procedure for odd / cases, which will be addressed in 
Section 4. Subsequently in Section 5 we apply our results to calculate thermodynamic quan- 
tities for dSp X S"^ and reproduce the first law of de Sitter thermodynamics. We conclude 
the paper in Section 6 and present numerical results for the warped branches in g = 5, 6 
and 7 in Appendix A. In Appendix B we describe the procedure to derive the perturbation 
equation used in Section 3. 

■^It is important to note that Kinoshita's ansatz corresponds to static sohitions, thus the perturbative 
solutions only indicates the existence of a static threshold mode. For the de Sitter solutions, this coincides 
with the onset of dynamical instability. However for the anti-de Sitter case, in order to have a physical 
dynamical instability, the perturbative modes must also violate the well-known Breitenlohner-Freedman 
bound; such considerations are beyond the scope of this paper. 
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2. Preund-Rubin-type compactifications 



2.1. The action 

We begin by reviewing the theory of Freund-Rubin-type compactifications. They de- 
scribe higher dimensional theories where the extra compact dimensions are stabihzed by the 
presence of a g-form flux. We will also include the presence of a bulk cosmological constant 
A. It was first observed by Freund and Rubin [1] that the presence of can induce a natu- 
ral compactification such that the solutions can be written as a product Aip x^/q, where Aip 
is the external p-dimensional spacetime with Lorentzian signature and A/'q is the compact 
^-dimensional space with Euclidean signature. 

The starting point for the theory is the following action (in units where G = c = 1) 

1 



dPxd'^y^{R-2A--FL , (2.1) 



where R is the scalar curvature, is the g-form field, and A is the positive cosmological 
constant. We will use lowercase Greek indices (/i, z/, . . .) to denote coordinates in the p- 
dimensional external spacetime, and lowercase Latin indices (m, n, . . .) to denote coordinates 
in the g-dimensional compact space. Uppercase Latin indices {M,N,...) will denote all 
coordinates from the entire (p + q) dimensions. 

Varying the action (2.1) gives the Einstein-Maxwell equations 

Rmn = 7 ■7rTRML2...LRN^^"'^'' 7^—.F'^9mn H ; -^Qmn , 

[q — ly. p-\-q — 2q\ p + q — 2 

In this work we are interested in solutions where the external spacetime is p-dimensional 
(anti)-de Sitter space, with a g-dimensional internal space. Following Kinoshita we consider 
the case where the internal space has SO{q) isometry. Such configurations can be described 
by Kinoshita's ansatz [8] 

(1 - h^r^) df+{l- h\^) dr^ + rMnf 



(P"2) 



+ e-^^*(^) [d9' + a{9r dnl_,^] , (2.3) 



where h is the Hubble parameter and df2^^^ denotes the metric of a fc- dimensional unit sphere. 
The g-form field that solves the Maxwell equations can be parametrized by its strength b and 
is given by = ha'^ e 9-2 a dr2(q_i). The warp factor $(6') and a{d) are functions 
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that depend only on one (finite) internal coordinate 9, which lies in the interval 9_ < < 6^. 
Loosely speaking, we may call 9- and 9+ as the 'south' and 'north' poles respectively. We 
further let 6'* be the position of the 'equator', whose definition will be made more precise 
below. In this case, the Maxwell equations are trivially satisfied, while the Einstein equations 
reduce to 

$" + (g - 1)— $' = (|)- l)/i2e *H \?Q 9-2 * e 9-2*, (2.4) 

a p + q — 2 p + q — 2 

„vi^ = -fi<f±l-H)a*-, (2.5) 

a [q — 2Y 

with a constraint equation 

/ ^^/ c^\f(^'^-^\ P(P + ?-2) ,2 / „ _2(p+g-2)^ 

(g-l)(g-2) ( = (g_2) ^ +P(p-l)^ e 

+ fe^e-^* _ 2Ae-^* , (2.6) 

where primes denote derivatives with respect to 9. The cosmological constant can be nor- 
malized to unity by rescaling b — 6A^/^, h — )■ hh}/"^, a — )■ aA~^/^ and 9 — )■ 6'A"-^/^; this 
rescaling implies that all length scales are measured in units of (anti)-de Sitter curvature. 

In the metric (2.3), there is still some freedom which is yet to be fixed, namely, the 
following transformations parametrized by a constant $0 

$^$ + $0, a^ae^*°, /i ^ /ie*% 6 ^ te^*", 

e^^e^*% t^te-*% (2.7) 

which leaves the metric invariant. In practice, this translates into a freedom in fixing the 
boundary condition for $ when integrating (2.4). Obviously this has no physical significance 
as the transformation in (2.7) leaves the metric and the physical observables unchanged. 
However, when performing numerical calculations, this freedom has to be exploited correctly 
to ensure that we discover all branches of solutions. In particular, a suitable choice of $0 is 
important so that the numerical analysis is able to discover the odd / branches (see Appendix 
A). 

2.2. The FR branch and its stability 

The Einstein-Maxwell equations (2.4) and (2.5) admit a trivial solution of the form 
(A)dSp X S*^, which we will refer to as the Freund-Rubin (FR) branch. In this case, the 
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solutions are 

Q 

$ = 0, am=pocos — . (2.8) 

Po 

Here, the north and south poles are 9± = ±po7r/2, and the equator is = 0. The components 
of the g-form field reduce to 

^^■■■^^ \ otherwise. ^ ' 

This solution is parametrized by three quantities: the Hubble parameter /i^, the fiux strength 
6^, and compactification radius po- We further note that the external space is de Sitter if 
h'^ > 0, Minkowski if h"^ = 0, or anti-de Sitter if h"^ < 0. The parameters {6^,/i^,po} are 
constrained by the Einstein-Maxwell equations. To see this we simply substitute the trivial 
solution (2.8) into Eqs. (2.4) and (2.5). This results in the following algebraic relations 

{p-l){p + q- 2)h^ + (g - 1)6^ = 2A , (2.10) 
(g - l)(p + g - 2)po 2 _ (p _ 1)6^ = 2A , (2.11) 

leaving one free parameter to describe the Freund-Rubin solutions. We can therefore describe 
this family of FR solutions as curves in the (6^, h?) phase space, which, according to the linear 
relation (2.10), are simply straight lines. 

The FR branch contains two channels of instabilities [5, 4]; the first is due to the 
homogeneous mode (/ = 0), which occurs for 

,9 2A , ^ 

b < 7 TT r • 2.12 

-(p-l)(g-l) ^ ' 

This instability is due to homogeneous perturbation of the compact sphere, which results in 
a uniform expansion or contraction of the compact space. Therefore, this instability does 
not change the shape of the compactification. 

The second instability, which is the main concern of this paper, are the inhomogeneous 
perturbations (/ = 2, 3, . . .) occuring when [4, 5] 

2A[/(Z + g-l)-2g + 2] 
- (p-l)[2(g-l)2-/(/ + g-l)]' ^ • ^ 

where — /(/ + g — 1) are the eigenvalues of the Laplace operator on S'^, corresponding to the 
eigenmode /. If the fiux satisfies (2.13), the linearized analysis show that there will be modes 
that deform the compact sphere which grow exponentially in time. The threshold stationary 
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g = 4 q = 5 q = Q 1 = 7 
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Table 1: The points of marginal instability, h'^/A) for p = 4 and various / and q 

mode is the limiting case where equality occurs in (2.13). This signals the presence of a new 
branch of solutions with deformed compact spheres according to the shape of its mode /. 
These points are called marginal stability points, and will be denoted by (6^, h^J. 

Thus, to obtain stable FR solutions, one must choose the parameters and di- 

mensions p and q to avoid the inequalities in Eqs. (2.12) and (2.13). For dimension q > 5, 
the inequalities (2.12) and (2.13) overlap for all h"^ and 6^, thus there are no more stable FR 
solutions for g > 5. 

From here onwards, we consider the inhomogeneous instabilities and the resulting new 
branches. The (6^, h"^) phase space gives us an intuitive interpretation of the stability of 
the solutions: when exceeds a particular threshold point indicated by Eq. (2.13), the FR 
solution becomes unstable to inhomogeneous perturbations. This leads to the possibility 
of a new class of stable solutions which are distinct from the FR branch. The intersection 
between the FR class and the new class should occur at the point of marginal stability 
(6^,/i^), where instability begins to arise. Table 1, based on the calculations of [1], shows 
various points of marginal instability for different cases of q. In particular, we can easily see 
that with the increasing number of compact dimensions q, the number of unstable modes 
increases as well.^ 

2.3. The warped branch 

For the case g = 4, Kinoshita has solved Eqs. (2.4) and (2.5) numerically and found 
that the warped branch is also a single-parameter family which intersects the FR branch at 
the point of / = 2 marginal instability. For convenience, we reproduce the phase diagram in 
Figure 1. We also produce phase diagrams for higher q cases in Appendix A for g = 5, 6 

^The points with negative values of are excluded as being unphysical. 
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q=4 




Figure 1: 6^ vs h"^ phase diagram of solutions q = 4:, first obtained by Kinosliita in Ref. [8] 

and 7 (Figure 2). Tlie boundary conditions required to solve Eq. (2.4) and (2.5) are: 

a(^±) = 0, \a'{9^)\ = l, $'(e±) = 0. (2.14) 

The zero of the function $ can be arbitrarily chosen due to freedom given by Eq. (2.7). In 
obtaining the solutions in Figure 1, Kinoshita made the following choice [S]: 

$(^*) = 0. (2.15) 

For the case of warped branches, we define the equator 6* = 6'* as the location where a'(^*) = 
0. Correspondingly a(6'^,) = p will be called the 'compactification radius', in analogy to the 
FR branch. 

3. Perturbative expansion of the Kinoshita ansatz: even / cases 

The warped branches obtained by Kinoshita in Figure 1, and the higher / modes shown 
in Figure 2 appear to be straight; thus it was conjectured in [8] that this branch may also be 
linear in /i^ and b"^, similar to Eq. (2.10). However, as we will see below, this is not the case 
as this branch is actually slightly curved. We will show this by obtaining an approximate 
analytical expression for the warped branch and show that there is a nonzero quadratic 
dependence on Sb"^, where 56^ is the amount of deviation away from the FR branch (this will 
be defined more precisely below). We first focus on the cases where / is even, i.e., / = 2N 
where is a positive integer. The perturbations in the odd modes (/ = 2N + 1) introduce a 
certain subtlety which we will address in Section 4. We will also derive analytical solutions 
for $ and a under the approximation where $ is small. This method reproduces the points 
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of marginal stability initially obtained by Bousso et al. [5] and Martin [1], therefore we shall 
refrain from assuming their results. 

In this treatment, the warped branches in (6^, h'^) phase space can be written as an 
expansion from a (yet unspecified) point (6q, h^) on the FR branch using 6b^ as the expansion 
parameter. Therefore we write 

b^ = bl + 6b\ = hl + wi6b^ + W2{6b^Y + ... (3.1) 

where the quantities b^, h^, wi, W2 ■ ■ ■ will be determined from the perturbed equations. 
Before proceeding to the details, we briefly outline the procedure that will be used to calculate 
wi and W2- 

• Firstly, we write $ and a as a perturbation series expanded about the FR branch (2.8). 

• In the first-order perturbation of (2.4), the requirement that the solutions are regular 
determines the points of marginal stability (6^, h"^). 

• The second-order perturbation of (2.4) results in an inhomogeneous differential equa- 
tions whose source term is parametrized by wi. Demanding the solutions to be regular 
constrains Wi to a particular value, giving the first order coefficient to the phase space 
equation (3.1). 

• The same step applied to the third-order equations will correspondingly give the value 
of W2. 

Our starting point is to expand Eqs. (2.4) and (2.5) about the FR solution (2.8). In the 
warped solutions, the compact space is no longer a perfect sphere, and therefore the com- 
pactification radius is no longer a constant. However, recalling the previous section, we 
define a quantity called 'compactification radius' as a(^*) = p, and thus the perturbation 
solutions will introduce higher order corrections in the form p = po + piSb'^ + p2 (56^)^ + . . . 

The problem we are considering in perturbation theory is to investigate the behavior of 
$ and a as the parameters (6^, h"^) move away from the FR branch in accordance to Eq. (3.1). 
We begin by writing $ and a as a perturbation series under the parameter 56^ about the FR 
solution: 

$ = + <t>idb^ + $2 [db^Y + $3 {Sb^Y + ... (3.2) 

a = Po cos h aiSb'^ + a2 [5b^Y + • • • (3.3) 

Po 
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Some care must be taken when choosing the appropriate coordinate system to use in the 
perturbation analysis. It is important to note that the coordinate 9 that appears in Eq. (3.3) 
has dimension of length, and hence will be physically deformed when the FR branch is 
perturbed. In particular, we may expect the locations of the 'poles' to be shifted by the 
perturbations. This introduces an ambiguity in terms of the coordinate locations where the 
boundary conditions are to be imposed. We address this in more detail in Appendix B where 
we derive the perturbation equations by introducing a suitable dimensionless coordinate 
system that allows us to impose boundary conditions at the correct locations. In the following 
subsections (3.1-3.3) we solve the perturbation equations from first to third order. The reader 
interested in the results may skip the details and proceed to Section 3.4 for the perturbative 
analytical solutions in (6^, h"^) space describing the warped branch. 

3. 1 . First- order perturbations 

The first-order perturbations of (2.4) is (see Appendix B) 



where we are using the coordinate x = sin — , and 



$1 = pIF, , (3.4) 



q-2 



Fi = (p - 1)^1 + ^ \ . (3.5) 
p + q - 2 

where we recognize Eq. (3.4) as an inhomogeneous Gegenbauer differential equation with a 
constant source term (see, e.g., [20, 21, 22]). Regular solutions exist if po and (3 satisfy the 
relation 

Pq/3 = /(/ + g — 1) , where / is an integer. (3.6) 

When the above equation is solved for 6q, and using Eq. (2.10) to eliminate h^, we obtain 

,2 ^ 2A [/(/ + g-l)-2g + 2] 

(p-l)[2g2_4g + 2-/(/ + g-l)] ' ^ ' ' 

which is precisely the saturation point of Eq. (2.13), i.e., the points of marginal stability 
6^. Therefore we have obtained the first result from our perturbation analysis: the warped 
branch, as specified approximately by Eq. (3.1), emanates from the FR branch at bo = b,^. 
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The explicit solution to Eq. (3.4) is given by 

$i(x) = A;iC/")(x) + ^, (3.8) 

where A = (g - l)/2 and c/^^ (x) is the Gegenbauer polynomial of degree / and order A. ki 
is an integration constant determined by the boundary conditions $(0) = 0. 

To complete the first-order analysis, we also seek the first-order solution of (2.5), which 
will be important in the subsequent sections that follow. To first order, (2.5) is 

(l-x^)^ + x^-ai = 0. (3.9) 
ax'' ax 

We see that (3.9) requires ai{x) to be even about the equator. This imposes a boundary 
condition on its first derivative, namely a'i{0) = 0. Together with ai(0) = pi, the solution 
that satisfies the boundary condition is 



1 (^V 1 — x"^ + X arcsin (x) j . (3.10) 



ai[x) = p 

The explicit value of pi can calculated by evaluating Eq. (2.6) at a{x)\x=o up to first order. 
Finally, using the above expression of ai(x), we can obtain the function which we call Xi in 
Appendix B (Eq. (B.2)), 



X, = -Pi arcsin {x)Vl^ = _lIEZ^ . (g.n) 
Po Po dx 



3.2. Second- order perturbations 

The value of Wi is determined by the requirement that the solution $2 of the second 
order perturbation must be regular at the poles. To obtain this solution, we expand (2.4) to 
second order by substituting Eq. (B.3) into (2.4) and collecting the second-order terms. 



X 



dx^ 



U 2 n 



$2 + ^1 



d$i 
dx 



2pg 



(g-2)2 
q-2 



{p~l)ip + q-2)hl + 



p^iq-1] 



p + q - 2 



2p^A 
p + q - 



p{q - 1] 



p + q 



Po 



X 



(p- l)(p + 9 

arcsin (x)\ d$i 



2)wi 



$1 + (p - l)PoU;2 



x^ 



dx 



(3.12) 
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Using the results obtained in the previous section for $i and Xi (Eq. (3.8) and (3.11) 
respectively), the second term in the parentheses on the LHS can be differentiated explicitly, 
and the resulting expression will cancel some terms on the RHS to give the following resulting 
equation 

d2 d " 



$2 = p^F2(x) + pi ( 2^Fi + (p - 1)W2 

Po 



(3.13) 



where 



Fo(x) 



(g-2)2 
2 



q-2 



{p-l){p + q-2)hl 
{p- l){p + q-2)wi 



p\q~l)^ ^2 
p + q-2 

piq-lf 



2p2A 



p + q - 2 



$2 



p + q-2 



Po 



(3.14) 



This is again an inhomogeneous Gegenbauer equation, where the source term F2{x) is now 
a function of x. This equation can be solved using the Green's function 



._v^ m+mxr (x)cf (xO(i-x-^)^-v^ 

'^"^'"^^ ^7r2i-2Ar(j + 2A) l{l + q-l)-j{j + q-l) 
which allows us to construct a complete solution to Eq. (3.13) as follows: 



(3.15) 



<l>2{x) = k2Cl^\x) + 



2piFi , {p-l)w2 



+ pI I dx' Giix,x')F2ix'). 



(3.16) 



Po /3 13 

Regular solutions exist only if F2{x) satisfies (see, e.g., [23] or any standard mathematics 
textbook) 



J dx F2{x)Ci^\x) {1 - x^) 



2XA-1/2 







(3.17) 



The remaining unspecified parameter in F2{x) is wi, which we may adjust to satisfy Eq. (3.17). 
For / = 2N, N = 1,2, ... we calculate Wi explicitly by direct integration of Eq. (3.17) using 
the various properties of Gegenbauer polynomials [20]. This gives 

hi + {q-l)Q 



Wi 



pbl-{p-l)ip + q-2)Q 



(3.18) 



where 

Q 



X 



-1)^ r(A + N)^T{2X + 3N)T{2N + 1)(A + 2N) 

4 r(iv + i)2r(A)r(A + 3iv + i)r(2A + 2Ar) 

[bl + {p- l)hl] [p\{q - l)hl + {p-l){p + q-2){2p + q- 2)H 



K[q-2Y{p-iy 



(3.19) 
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Before concluding this section, we complete the second-order analysis by obtaining the 
solution of 02- This solution will be required in the following section. Equation (2.5) to 
second order is 



da. 



X 



X- 



0-2 



dx^ ' dx 

pI / (arcsin (x))^ 

Po V Vl - x^ 



VT 



X' 



(3.20) 



Unlike the earlier case of ai, the equation will be different for each particular /, because the 
source terms depend on $i. The second-order equation is solved by 



a2(x) 



CiX + C2 



— x^ + ^arcsin (x) 



TT 



X 



arcsm (xj v 1 — x^— — 

2po dx 



p(p + g-2) 
(.-2)^ 



Vl — x^ + X arcsin (x 



dx X 



X arcsm x , 



d<l>i 

dx 

d$i 

dx 



(3.21) 



Here c\ and C2 are the integration constants to be determined from the boundary conditions. 
Similar to the first-order case, Eq. (3.20) is symmetric about the equator, which gives the 
same boundary condition 02(0) — 0- The second integration constant can be fixed using 
02(0) = P2- 

Finally, the function X2 will accordingly be the second-order term of the expansion of 
Eq. (B.l). The specific form will depend on the particular branch, because it depends on 
ai, 02, and $1; therefore will be too cumbersome to be shown here. Nevertheless it is a 
straightforward calculation once the expressions for ai, 02, and $1 are at hand. 



3.3. Third- order perturbations 

It is evident by now that to analytically obtain the coefficient Wi, it is necessary to solve 
the perturbation equations up to {i + l)-th order. To obtain the third-order equations, we 
proceed in a manner completely analogous to the procedure in the second-order perturba- 
tions. We substitute (B.3) into (2.4) and collect the terms which are cubic in 6b^ to produce 
the third-order analogue of Eq. (3.12). This equation will depend on the functions $1, $2? 
ai, 02, Xi, and X2, which have been obtained in the previous sections. Substituting these 
results into the third-order equations gives the corresponding analogue of Eq. (3.13), after 
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some simplification, 

^ '^^^2" ~ " Po^sW - Po ( ^1 + —yP- ^)W2 ) , (3.22) 

where 



d$i 

da; 

2f{x)VT^\ ,^ , , + 9 - 2) (1 - ^d$i 

— P<I>1 



F3(x) = ^ /^^^ _ ^ ^ ^ /^^^ ^ ^(arcsin {x))^Y^ 

6 V Po / V Po 



+ 



pi^x^ ) (g — 2)2 PqX \ dx 

/2(p^l)(p + ,-2)^,^ + !«B3 + H««±ii;3) 
V ?-2 Po Po / 

1^ ,P1rA^2 2/(x)v/r^ d$i 

-54 + $1 ^ , (3.23) 



with the following constants abbreviated as 



51 = 


(p- 






(P- 


^)( 


^3 = 


(P- 


^)( 


^4 = 


(P- 





2(p + g- 


2) 


g-2 




2{p + q- 


2) 


q-2 




2{p + q- 


2) 


q-2 




2{p + q- 


2) 


q-2 



p + g-2V g-2 y p + g-2Vg-2^ 

' ,-1 ppto-i)V,._ ---^^ 



p+g-2V g-2 y " p+g-2Vg-2 
g — 1 / 2p{q — 1) 



p + q — 2 \ q — 2 

\,,^(Mi^V. (3.24) 

p + q-2 \ q-2 



PoVT^ J dxx (^^^ • (3.25) 



and the function f{x) is given by 

_ p{p + q- 2) 
(g-2)2 

The particular solution $3 can be evaluated analytically by a Green's function integral similar 
to Eq. (3.17), where for regular solutions, we require 

j' dx F3{x)Cl^\x) (1 - = , (3.26) 

giving the value of W2- 



3.4- Results 

Summarizing the results on the previous sections, we find that the phase space of the 
warped branches can be approximately described by 



(3.27) 
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q: 


4 


5 


6 


7 


I = 2: 


7 

54 


431 
3444 


199 
1668 


8 
71 


I = 4: 






2293 
10056 


646 
2547 



Table 2: Exact values of wi for p = 4 and various g, for even / branches 





q: 


4 


5 


6 


7 


/ 


= 2: 


722 


315975 


91208835 


197372 


216513 


9790144 


1117217504 


1073733 


/ 


= 4: 






24624514905 


100489280 






3041267079296 


1835880147 



Table 3: Exact values of W2 for p = 4 and various g, for even / branches 

where the values wi and W2 are derived in the previous sections. We list the results in Tables 
2 and 3 for p = 4, and the cases g = 4, 5, 6 and 7. 

Note that the values of W2 are nonzero, which shows that the warped branch is not 
purely linear in the {b'^,h'^) phase space. For example, in the case studied by Kinoshita 
{p = q = A, I = 2) , the warped branch is (in units where A = 1) 

For this particular case, we have /i^ = 1/18 in the first term which is in expected agreement 
with [5, 4]. The next two terms are with wi = —7/54 and W2 = —722/216513 which are 
precisely obtained from Eqs. (3.18) and (3.26) respectively. It has been checked that the 
results of Tables 2 and 3 agree with the numerical results within numerical accuracy. 

4. Odd / cases 

For the odd / cases, if we apply the procedure of the previous section, it may appear 
that the perturbative solutions to Eqs. (2.4) and (2.5) could not be determined uniquely. 
The reason for this can be seen from Eq. (3.8). If / is odd, then Ci^\x) is anti-symmetric, 
thus vanishing identically at x = 0. This implies that at the equator, $i(x = 0) = Fi/ (3^0, 
which is inconsistent with the choice in Eq. (2.15). 

To overcome this problem, we exploit the freedom given by Eq. (2.7) and shift the 
zero point of $ by choosing $(x = —1) = (this is equivalent to $(^-) = 0, i.e., $ is 
zero at the south pole) and perform the perturbation analysis in this scheme. The results 
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Q- 


5 


6 


7 


/ 
I 


= 3: 
= 5: 


174231551 
577251843 


56669459 
166268616 


5316380 
15099489 

5029755010667 
24287636988414 



Table 4: Exact values of wi for p = 4 and various g, for odd / branches 



can be consistently transformed back via the inverse of Eq. (2.7) to obtain the approximate 
solutions for the warped branches. Such transformations must be handled with care; the 
transformation (2.7) involves a rescaling of b^. Since 6^ is our perturbation parameter, 
different perturbation orders will be mixed up by the transformation (2.7). 

Indeed, we find that information from the second-order perturbation results under the 
scheme $(^-) = is transformed into first-order under the scheme $(6'*) = 0. Therefore 
perturbations up to third-order, if converted back to the choice $i(x = 0) = Fi/P ^ only 
gives values up to Wi. To obtain W2 for odd / cases one may need to go beyond third-order 
perturbations, which is presently beyond the scope of our analysis. The numerical values for 
Wi in the odd / branches are given in Table 4. For example, in the case p = A, q = 5 the 
phase diagram for the warped branch is given by (in units where A = 1) 

22 174231551 



h' = 6b' + O 

315 577251843 



W • (4.1) 



Here we have /i^ = 22/315, in agreement with [5, 4], and the next term is with wi 
-174231551/577251843. 



5. Thermodynamics 

As an application and consistency check of our perturbation results, we turn to the 
study of de Sitter thermodyamics. Below we will use our perturbation equations to derive 
the first law of thermodynamics, which was first derived by [18] using a different method. 
Clearly, the following analysis will only apply to solutions containing a horizon, i.e., h"^ > 0. 

In [18], Kinoshita and Mukohyama defined the entropy S and total flux J-" as the basic 
thermodynamic variables: 

S = iV^)%i) / d^ a^- , (5.1) 
J^= (f = / d^ a^-ig-^^* . (5.2) 
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By a scaling argument and variation of the action (2.1) (see Ref. [18] for details), the following 
first law of thermodynamics was derived: 



dS 



n 



{p-2)t 



-dJ-. 



(5.3) 



A{p - l)hP 

We can show that Eq. (5.3) can also be derived from the perturbation analysis. We consider 
the changes in S and J-" as the warped solution moves away from the marginal stability 
points. Hence the parameters changes according to Eq. (3.1), which we rewrite here for 
convenience: 



6b' 



+ U)i5b'' 



(5.4) 



Here the subscripts in i5w and Jv indicates that we are referring to the entropy and flux for 
the warped branch. Recall that when 56^ = 0, the warped and FR branches coincide, so at 
Sb'^ = 0, we have 5w = i5fr = Sm and = J-'fr = J^m- For 56^ ^ 0, small increments to 
Sb'^ induces corresponding changes to the entropy and flux. 



(5.5) 



To calculate 5S^ and 5J^^, we expand (5.1) and (5.2) to first order and make use of the 
results (3.8) and (3.10). Evaluating the resulting integrals, we find, to first order 



5S^ 



^(p-2)^{g-i) v^r (I) 



.2/ J. ^^Po 2(p + g-2)Fi {p - 2)wr 
,(^TyPo|(.-l)- ^ 



+ (g-l 



Pi v^r(g-l) 

Po2-ir(|)r(^) L 



6b' 



(5.6) 



where ip{x) is the Euler psi function defined by ip{x) 
order change in flux is 



dx 



Inr(x) [20]. Similarly, the first 



6m^(g-l) 



r(^) 



2bl 



po 2p(g-l)Fi 
Pi 



+ 



Pi v^r (g - 1) 

Po2^/-ir(|)r(^) 



q+ 1 



/3 
q-l 



6b' 



(5.7) 



Eliminating 6b'' (with the aid of Eq. (2.10)) Eqs. (5.6) and (5.7) reduces to (5.3), reproducing 
the first law (up to first order, we have 6S = dS and 6J^ = dj-"). This demonstrates the 
consistency of the perturbative analysis with the first law of thermodynamics. Due to the 
fact that both the FR and warped branches obey Eq. (5.3), it follows from Eq. (5.5) that at 
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first order, the FR and warped solution with the same flux also have the same entropy. It 
is reminiscent of the analysis of black brane stability, where for the same mass, the uniform 
and non-uniform black branes have the same entropy up to first order [19]. 

If we extend the comparison by including second-order terms, we find that 



, (5. 



which means the respective entropies at second order are still the same. When the difference 
in the above equation is calculated, there is a nonzero term of order (56^)^, which may 
possibly indicate that the entropy values for the different branches begin to differ at third 
order. However to verify this one has to calculate the expansions of S and J-" up to third 
order, which requires the exact solution for a^{x), which we do not have. 

We may compare the above analytical calculations with the numerical results of [IS]. 
The entropy difference iSw — iSfr for fully non-perturbative (numerical) solutions of Eqs. (2.4) 
and (2.5) was calculated in [18] and shown to be consistent with the correlated stability 
conjecture [17], i.e., iSw > iSfr when the FR branch is dynamically unstable. An inspection 
of Figure 3 of Ref. [18] shows that iSw — iSfr as a function appears fiat for regions close to 
the marginal stability points. Our calculations above identify the scope of this region to be 
at least second order in deviations away from the marginal stability point. 



6. Conclusion 

In this work, we have used Kinoshita's ansatz to find numerical solutions of various 
modes for the cases g = 4, 5, 6 and 7. It was shown that the warped branches under this 
ansatz are consistent with the stability analysis for the FR branch by Martin [ 1] and Bousso 
[5] . In particular, the points of marginal stability gives rise to a new branch emanating from 
the FR branch. Although the warped branch in phase space appears to lie on a straight line, 
perturbation theory shows that this is not true because W2 is indeed small, but still nonzero. 
The warped branch has a quadratic dependence and is not purely linear in b^. 

Using perturbation theory, expanding Kinoshita's equations of motion about the FR 
branch reproduces the points of marginal stability (^mi^m)- Also we were able to derive 
the explicit solutions of $ and a up to second order. Furthermore it was found that the 
requirement that the solutions of $j must be regular provides constraints between h"^ and 
b^, which gives the approximate equations which describe the phase diagram for the warped 
branches in regions close to the FR branch. 
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As an application to de Sitter thermodynamics, we calculated the thermodynamic vari- 
ables S and J-" for the warped branches and reproduced the first law of thermodynamics. 
Our calculations show that for the same flux J-", the two branches have equal entropy up 
to second order. However it is already known from the numerical results of [IS] that the 
warped branch has higher entropy when the FR branch is dynamically unstable. Therefore 
we expect that the entropy difference between the two branches should enter at third-order 
or beyond. 

In addition to the above mentioned results, we have also found numerical evidence for 
the existence of yet another new type of branch which is not yet captured by any previous 
study. This branch intersects with the g = 4, / = 2 branch, but does not intersect with 
the FR branch. It has less symmetry than the / = 2 branch (it is neither symmetric nor 
anti-symmetric about the equator) and satisfies the boundary conditions defined above. The 
intersection with the warped q = 4, I = 2 branch occurs around the region h"^ ~ 0.0199A, 
~ 0.0951A, which is where the / = 2 is unstable and hence beyond the perturbation 
analysis of Ref. [IS]. Although not much is known about this branch at this point, it indicates 
the possibility that warped flux compactifications exhibit a rich structure and there may be 
many families of solutions which might deserve further study. 
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A. Numerical results for g = 5, 6 and 7 

Eqs. (2.4) and (2.5), in general, can be solved for any choice of p and q. We obtain the 
solutions for the cases p = 4, 4 < g < 7. For each case, the warped branch intersects with 
the FR branch at the point of marginal stability (6^, h"^). See Figure 1 for the g = 4 case 
[8] and Figure 2 for g = 5, 6, 7. 

In order to obtain all possible branches, the boundary conditions used to solve the 
Einstein equations are Eq. (2.14), which are slightly different from the one used by [S].^ In 
our case, the Einstein equations are integrated from the north to south pole, and hence, 

^In [(S], one of the boundaries is the equator, and the condition imposed is $'(6'*) = 0; which will produce 
the even I solutions, but not solutions with odd I. 
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there are no constraints imposed on $ at the equator, and thus solutions with $'(^*) = 
(even / solutions) and $'(6'*) 7^ (odd / solutions) can be found. This is a straightforward 
generalization of the numerics of [8] which allows us to find the odd / solutions which were 
previously not considered. 

The specific shape of the solutions $ will have similar shapes to the particular / modes 
from which it emanated from the relevant FR instability point, thus we name each branch 
by their corresponding /'s. As we have seen from Section 3, at regions close to the FR 
branch, the warped solutions $ are approximately Gegenbauer polynomials plus correction 
terms. Taking a specific example, in g = 6, the / = 3 branch is the one which emanated 
from the / = 3 instability 6^ = 14A/39, = A/117. Close to the FR branch, $ is pro- 
portional to 6*3^^^^ (sin ^). In the non-perturbative regime, $ will still inherit some features 
of C3^^^''(sin ^), namely it is anti-symmetric about 6^, and has four turning points. On the 
other hand, the solution for a in all cases have the same general appearance as shown in [S]. 

B. Establishing the perturbation equations and boundary conditions 

To derive the perturbation equations of (2.4) and (2.5), we need to establish a suitable 
coordinate system to remove the ambiguity of the boundary locations. Our starting point 
is to first consider the FR branch described by Eq. (2.8). It is intuitively clear that 0/pq 
is a dimensionless angular parameter that labels the points along one direction of the q- 
sphere with radius po- This also means that 6 has dimension of length and takes values 
-vrpo/2 <e< TTpo/2. 

Since we are looking for a warped solution by perturbing the FR branch, the length 
scale associated with 6 will change when the FR solution is deformed. This introduces an 
ambiguity in the endpoints of 6 (which are the location of the poles). This is problematic be- 
cause there is now an ambiguity in the location for which to impose the boundary conditions. 
To approach this perturbation systematically, it is convenient to define 



Since for both the warped and FR branches, boundary conditions require a{6) to vanish at 
the poles. So the poles can be specified precisely as a{X)\x=±i = 0. For the case of the 



a- = ctwarpcd is cxpectcd to be the corresponding generalization to include the warped branch. 
Recall that when perturbing the FR branch to find approximate warped solutions, a{6) is 




(B.l) 



FR branch, apii{0) = Pocos— , we simply have X = sin—. The coordinate X for the case 
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q=5 



q=6 




Figure 2: 6^ vs phase diagrams for cases g = 5 (top left), g = 6 (top right) and q = 7 
(bottom). The smaller branches for the cases q = 6 and q = 7 are magnified in the insets 
for clarity. 



written as a series in Eq. (3.3). Inserting Eq. (3.3) into (B.l) and expanding will give X as 
a series which we denote by 

X = po sin — + X^6b^ + X2 {6b'') ' + ..., (B.2) 

Po 

where Xi,X2, . . . are the coefficients of the respective orders of the expansion. We find it 
convenient to use x = sin — as our coordinate in the differential equations. With the inverse 
relation 6 = po arcsin(x), the X defined in Eq. (B.l), and its expansion (B.2) can be viewed 
as functions of x. Therefore the perturbation series of the warp factors can be calculated 
with X as a variable (—1 < x < 1), with corrections due to the shifting of the poles correctly 
taken into account when X is evaluated up to appropriate order. 
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Using our new coordinate system, we can now write the series in Eq. (3.2) as follows, 
$ (X) = $ (^x + + X2 {6b^Y + . . .) 

d$i(x) 



d$i(a;) 
dx 



+ $3(^)+^l 



d$2(x) 
da; 



dx 



2 ^ da;2 



{5b^Y + {5b'') 



(B.3) 



where the second equality follows from Taylor expansion up to third order in 6b'. Then, to 
derive the perturbation equations, we insert (B.3) into Eq. (2.4) and expand to appropriate 
order. The procedure is completely analogous for the perturbations of a. 
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